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ABSTRACT 



We study the gravitational lensing magnification produced by the intervening 
cosmological matter distribution, as deduced from three different hierarchical 
models (SCDM, LCDM, CHDM) on very high redshift sources, particularly 
supernovae in protogalactic (Pop III) objects. By means of ray-shooting 
numerical simulations we find that caustics are more intense and concentrated 
in SCDM models. The magnification probability function presents a moderate 
degree of evolution up to z 5 (CHDM) and z ^ 7 (SCDM/LCDM). All 
models predict that statistically large magnifications, /i ^ 20 are achievable, 
with a probability of the order of a fraction of percent, the SCDM model 
being the most efficient magnifier. All cosmologies predict that above z ^ A 
there is a 10% chance to get magnifications larger than 3. We have explored 
the observational perspectives for Pop III SNe detection with NGST when 
gravitational magnification is taken into account. We find that NGST should 
be able to detect and confirm spectroscopically Type II SNe up to a redshift of 
^ fa 4 in the J band (for Tsn = 25000 K); this limit could be increased up to 
z ~ 9 in the K band, allowing for a relatively moderate magnification. Possibly 
promising strategies to discriminate among cosmological models using their GL 
magnification predictions and very high- 2; SNe are sketched. Finally, we outline 
and discuss the limitations of our study. 



Subject headings: Cosmology: theory - dark matter - gravitational lensing - 
supernovae - galaxy evolution - methods: numerical - statistical 



-3- 



1. Introduction 

Current models of cosmic structure formation based on CDM scenarios predict that 
the first collapsed, luminous (hereafter Pop 111) objects should form at redshift z ^ 30 
and have a total mass M ft; lO^M© or baryonic mass ~ lO^M© (Couchman & Rees 
1986, Ciardi & Ferrara 1997, Haiman et al. 1997, Tegmark et al. 1997, Ferrara 1998). This 
conclusion is reached by requiring that the cooling time, tc, of the gas is shorter than the 
Hubble time, tn, at the formation epoch. In a plasma of primordial composition the only 
efficient coolant in the temperature range T < 10^ K, the typical virial temperature of 
Pop III dark matter halos, is represented by H2 molecules whose abundance increases from 
its initial post-recombination relic value to higher values during the halo collapse phase. 

Thus, as the collapse proceeds, the gas density increases and stars are hkely to 
be formed. However, the final product of such star formation activity is presently 
quite unknown. This uncertainty largely depends on our persisting ignorance on the 
fragmentation process and its relationship with the thermodynamical conditions of the 
gas, both at high z and, in a less severe manner, in present day galaxies. Ultimately, this 
prevents firm conclusions on the mass spectrum of the formed stars or their IMF. This 
problem was already clear more than two decades ago, as pioneering works (Silk 1977; 
Kashhnsky & Rees 1983; Palla, Salpeter & Stabler 1983; Carr et al. 1984) could not reach 
similar conclusions on the typical mass range of newly formed stars in the first protogalactic 
objects. Roughly speaking, two possibilities can be envisaged: either (i) a limited number 
of Very Massive Objects (VMOs), single stars with mass in the range 10^ — lO^M© could 
be formed, or (ii) a more common stellar cluster, slightly biased towards low-mass stars (or, 
even, some combination of the two involving low- mass star coalescence to form a VMO). 
Apart from the fact that observational evidences of the VMO hypothesis are still lacking in 
the local universe, there are also several theoretical problems (El Eid et al. 1983; Ober et 
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al. 1983) that make VMOs as unlikely stellar prototype candidates in the early universe. 
Then, if the IMF is more of the standard type, questions arise about its shape and median 
value. A common claim in the literature is that the absence of metals (as it is the case in a 
collapsing Pop III) should shift the peak and the median of the IMF toward higher masses. 
Obviously, assessing the presence of massive stars that produce ionizing photons and die as 
Type II supernovae (SNe) would be of primary importance to clarify the role of Pop Ills in 
the reionization and rehating of the universe, and, in general, for galaxy formation. Thus, 
studying the formation of the first stars has become one of the most challenging problems in 
physical cosmology. It is not the aim of this paper to investigate this aspect further; detailed 
studies can be found in Gnedin & Ostriker 1997, Haiman & Loeb 1997, and Ciardi & Ferrara 
1997. Here, we simply point out that if SNe are allowed to occur at these very high redshifts, 
they could outshine their host protogalaxy by orders of magnitude and likely become the 
most distant observable sources since the QSO rcdshift distribution has an apparent cutoff 
beyond z ^ A. Nevertheless, even for (future) large telescopes as NGST, VLT and Keck this 
will not be an easy task, this requiring to reach limiting magnitudes 33 in the near IR 
to observe a SN exploding at z = 10. Fortunately, the matter distribution in the universe 
behaves like a gravitational lens that distorces, and, most importantly here, magnifies the 
source flux. This magnification, as we will see, strongly enhances the detection probability 
of Pop III SNe, bringing their apparent magnitudes inside the instrumental capability range. 
A different, but complementary, strategy would be to monitor the central regions of nearby 
galaxy clusters (Miralda-Escude & Rees 1997), where magnifications high enough to allow 
for the detection of faint objects could be produced. We have chosen to investigate the 
gravitational lensing (GL) magnification effects of the entire cosmological mass distribution 
to verify if the different GL magnification predictions of the most attractive cosmological 
models could be used as a test bench of their descriptive power of the universe. Pushing 
this comparison to high redshift, as we will see, is crucial since magnification properties 
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deduced from hierarchical models start to differentiate considerably only beyond z 5. 

Gravitational lensing is sensitive to the cosmic mass distribution, which in turn 
depends on the adopted fluctuation spectrum, to the expansion factor, and to cosmological 
distances; therefore, it represents a perfect tool for cosmological studies. As stated above, 
our aim is to calculate GL effects on high z sources, as for example SNe in Pop III objects, 
and, in particular, the magnification probability as a function of the source redshift and 
of the cosmological model. In order to obtain this result it is necessary to define how- 
structure is formed in the universe and a method to evaluate the statistical incidence of GL 
effects. The former problem may be approached by the linear perturbation theory of cosmic 
density fluctuations, numerical N-body simulations or by semi-analytical methods, as, for 
example, the Press & Schechter (1974, hereafter PS) formahsm. GL can be investigated 
both by numerical simulations following light propagation or via semi-analytical estimates. 
We have chosen to solve the selected problem by using simulations based on the so-called 
ray-shooting numerical scheme as far as GL is concerned and to derive the mass and redshift 
distribution of the lens ensamble from PS. The reasons for this choice arc motivated below 
(see § 2); in general, we find that the ray-shooting -|- PS approach constitutes a fiexible and 
accurate enough tool for our purposes. The PS approach has been used also in conjuction 
with analytical techniques by several authors (Narayan & White 1988, Kochanek 1995, 
Nakamura & Suto 1997). 

This approach is somewhat different with respect to others used in past literature to 
investigate analogous problems. Schneider & Weiss (1988) have applied the ray-shooting 
scheme to study the propagation of light in a clumpy universe and to derive information on 
the infiuence of matter inhomogeneities on the apparent luminosity of lensed sources. The 
main difference with that work, is that we allow the lens mass to be distributed according 
to the prescriptions of a given cosmological model, rather than assuming identical clumps 
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of mass M. Essentially for the same reasons, our work differs also from the one presented 
by Lee & Paczynski (1990), who consider GL in a fiat spacetime and approximate the 
true three-dimensional distribution of matter with multiple identical screens. Jaroszyhski 
(1991) allows for density fluctuation to grow starting from the particular initial power 
spectrum oc (which is at the limit of the hierarchical structure formation model) and 
superimposing a galaxy population whose mass distribution is inferred from the Schechter 
function and a prescription for the M/ L ratio. The delicate assumption that galaxies trace 
mass in the universe, as well as the use of the Schechter distribution at high redshift, 
appears difficult to be justified, however. 

Massive computational methods coupling the full power of N-body simulations and 
ray-shooting start to become available in the literature (Wambsganss et al. 1996; Premadi 
et al. 1997), although with some unavoidable limitations which make them suitable for 
more focussed studies. Wambsganss et al. (1996) examined the evolution of various GL 
statistical indicators {i.e. multiple image frequency, lens redshift distribution, magnification 
maps, caustic and critical line geometry) in a CDM model up to z = 3. In the second 
work, gravitational shear and magnification of a single light beam produced by lenses in 
the redshift range 2; = — 5 for three cosmological models are calculated. Our work is 
intended to mainly explore the magnification effects produced by the cosmological matter 
distribution deduced from three different hierarchical models (SCDM, LCDM, CHDM 
whose parameters are given in § 2) on very high redshift sources. Particular emphasis 
is devoted to the perspectives offered by GL for the observational detection of SNe in 
protogalactic objects. We also explain how such observations could help constrain the 
validity of the proposed classes of cosmological models. 

The plan of the paper is as follows. In the next Section we define the cosmological 
models under study and derive their predicted lens mass distribution; in § 3 we introduce 
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and discuss the basic concepts of thick gravitational lenses. Sec. 4 is devoted to the 
description of the numerical scheme used for the GL simulations whose results are presented 
in § 5. Sec. 6 contains the observational implications of our results for future detection of 
high- 2; SNe with NGST; finally, the discussion and summary in § 7 conclude the paper. 

2. Cosmological Lens Distribution 

As stated above, our main aim is to determine the GL effects due to the intervening 
cosmological mass distribution on distant sources. Currently, the most attractive models 
for structure formation in the universe, assume that perturbations grow from an initial 
gaussian field of fluctuations with a scale-free power spectrum. These perturbations become 
nonlinear at redshift ^50 and form bound objects, which subsequently merge to form 
larger and larger structures in a hierarchical fashion. In spite of the relatively widespread 
agreement on this general scenario, a detailed comparison of the predictions of such type 
of models with observational data does not allow yet to discriminate among the various 
hierarchical models so far proposed and to uniquely fix their free parameters. 

Given these uncertainties, and with the hope to constrain at least some of the properties 
of the proposed models, we will consider and compare three different cosmological models 
with the same total density parameter fl = Qm + f^A + flu — 1, contributed by cold 
dark matter -|- baryons, cosmological constant A, and hot dark matter, respectively. In 
particular, (i) Standard Cold Dark Matter (SCDM), has Qm — — — 0; (ii) Lambda 
Cold Dark Matter (LCDM) has Qm = 0.4, = 0.6, fi^ = 0; finally (%%%) Cold+Hot Dark 
Matter (CHDM) parameters are = 0.7, = 0, = 0.3. 

In order to specify the spectrum completely we have to fix both the present value of the 
Hubble constant Hq — 100/i km s~^ Mpc~^, and the normalization of the spectrum, often 
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parameterized by its gaussian variance (see eq. Q below) erg = cr{8h~^ Mpc). The value 
of h remains uncertain by about a factor of two: 0.4 < h < 1. Standard relative distance 
methods using Type la SNe (Garnavich et al. 1997, Perlmutter et al. 1998, Kim et al. 1997) 
tend to favor higher values oih 0.6 — 0.8, whereas methods based on fundamental physics 
methods, including GL, suggest lower values h ^ 0.4 — 0.6 (Grogin & Narayan 1996; Kneib 
1998 and references therein). For this reason, we adopt the intermediate value h = 0.65 
throughout the paper; fixing h has the advantage that, when comparing the predictions of 
the various models, one can isolate differences that depend only on their intrisic properties 
listed above. The second point concerns the normalization of the power spectrum. The 
actual form of the power spectrum can be obtained by numerically integrating the relevant 
numerical equations; we use the analytical fit to such results provided by Efstathiou et 
al. (1992): 

I e |2 /-I \ 

I '^l ^ 7 r 77o „ii/. 2/1/ ' 

{l+[Bk + iCkY/^ + iDkf] } 

where i/=1.13; 5=6.4 hT Mpc; C=3.0 hV Mpc; D=1.7 hV Mpc, n=l; T is the shape 

parameter and it is equal to T = h for SCDM, T = hVLu for LCDM, V = 0.2(0.3/r],,)^/2 

CHDM. The coefficient A, or equivalently the value of as, must be determined by choosing 

an appropriate normalization. We fix ag such that the number of clusters present at 2; = 0, 

about one per {lOOh^^MpcY , is correctly reproduced. This implies = 0.64, 1.07,0.68, for 

SCDM, LCDM, CHDM, respectively. 

Once the spectrum has been fully specified, the structure formation and evolution 
can be determined, for example, via N-body numerical simulations. Here we take a 
complementary approach that makes use of the Press-Schechter (PS) formalism (Press 
& Schechter 1974; Bardeen et al. 1986; Peacock & Heavens 1990; Bond et al. 1991) in 
its standard form; we will briefiy discuss later in this Section the main advantages and 
limitations of such approach. Given a power spectrum (in our case eq. H) one can 
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write the gaussian variance of the fluctuations on the mass scale M: 



^li = J^,W\k,R)\5,\\ (2) 

where M = (4/3)7rpi?^, p is the matter density, and 

W = [smikR) - (kR) cos(fci?)] ; (3) 

is a top-hat fllter function. From the results of the nonlinear theory of gravitational collapse, 
stating that a spherical perturbation with overdensity 6c = Sp/p > 1.69 with respect to the 
background matter collapses to form a bound object, through the PS formalism we can 
derive the normalized fraction of collapsed objects per unit mass at a given redshift: 



/(M..) = ,/fMi+i),-.?(™,(_^|. (4) 



A useful quantity is the cumulative mass function 



F{M, z) = I f{M, z)dM = erf ^ ^'^'^J j . (5) 







In Fig. 1^, as an example, we show the function f{M,z)M, which gives the fraction of 
collapsed objects per unit logarithmic mass for the three cosmological models and two 
different redshifts. 

The PS formalism has been shown to be in surprisingly good agreement with the results 
from N-body numerical simulations (Brainerd & Villumsen 1992, Gelb & Bertschinger 1994, 
Lacey & Cole 1994, Klypin et al. 1995). Nevertheless, it suffers from several limitations 
(for a discussion see, for example, Yano et al. 1996; Porciani et al. 1996), and that the 
complexity of structures in the universe cannot be described in full detail. There are at 
least two pieces of information that are not present in the PS recipe: (i) the internal 
structure, i.e. the density proflle, of the collapsed objects, and (ii) their spatial distribution 
at a given redshift. In general, these aspects are relevant in GL experiments and for this 
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reason, as we will discuss in the next Section, one has to make some additional hypothesis. 
Since we approach the study of GL effects using statistical methods, the bias introduced 
by those assumptions needs to be evaluated carefully; we leave this to the final Section. 
However, obtaining the mass distribution as a function of redshift with PS, not mentioning 
its widespread use in literature and in GL studies that allows for straight comparisons 
and checks, has the advantage that it is computationally convenient and therefore the 
parameter space exploration and the optimization of the GL numerical scheme (see §4) 
can be performed with relatively unexpensive procedures. For these reasons we feel that is 
worthwhile to attack the problem with these approximate techniques that can represent a 
useful complement to more massive future computations. 



3. Thick Gravitational Lenses 

In this paper, we are interested in determining the GL magnification properties of high 
redshift sources (typically SNe) due to the intervening cosmological matter distribution. 
Since this matter distribution, as obtained from N-body simulations and/or PS formalism, 
covers a large range of redshifts, the standard assumption that the lensing mass is 
concentrated on a single plane perpendicular to the propagation of the reference light ray 
{i.e. the optical axis) has to be modified to include such spread. This can be achieved by 
splitting the mass distribution on a series of planes, and furthermore postulating that (i) 
deflections on the various planes are independent and (ii) the trajectory of a light ray is 
modified only at planes. These deflections can be calculated through the deflection potential 
"^zi^z), where z refers to the redshift of the plane and is the impact parameter on that 
plane; the propagation between planes can be generally described by the Dyer-Roeder 
(1972) equation, also discussed by Futamase & Sasaki (1989), Babul & Lee (1991), and 
Nakamura (1997). In brief, the main requirement when using Dyer-Roeder distances is 
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that the redshift intervals are ^ 1, in order to fulfill the underlying hypotheses of the 
theory, i.e. a light ray bundle is not affected by nonlocal gravitational effects (shear) and 
significative local deflections. The validity of the above scheme, which pundits call a "Thick 
Gravitational Lens", is discussed in detail by Seitz et al. (1994), Seitz & Schneider (1994), 
but also in Schneider et al. (1992) and Schneider (1997). In the following we describe how 
we model the lens systems corresponding to the various cosmological models using the 
above approximation. 

Let us suppose to split the mass on planes, at redshift ,21,-22, ..,^n5 --j^at, equally 
spaced with redshift interval A2; < 1 both for the above reasons related to the use of 
Dyer-Roeder equation, and to correctly sample the cosmic mass distribution eq. ^. A sketch 
of such conflguration is given in Fig. ^ Then the impact parameter on the plane n + 1 is 
given by 

t (1 + ^n-l)-Dn,n+l ^ -0^-1, n+1 ^ ^ y n /^n 

U+1 — TT- H U - ^n,n+iy S.^ zA^n), \P) 

{i- + Zn)jJn-l,n J^n-l,n 

where the deflection potential is 

S2 is the mass surface density on the plane at redshift z and is an arbitrary length scale. 
The angular distance Dij = D{zi,Zj) is determined by solving the Dyer-Roeder equation in 
the form suggested by hinder (1988a, 1988b) 

D + ^j^ID + + z)-'D ^(1 + s)As{z)Qs{z) = 0. (8) 

Here the dot indicates redshift derivatives; s is deflned by the equation of state p = E^sp^ 
involving the total pressure and the various density contributions, including the cosmological 
constant, in a Friedmann-like universe; for example, s = — 1 for the cosmological constant, 
whereas s = for a dust universe. Finally, a convenient expression for the generalized 
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deceleration parameter q{z) is (Linder 1988a) 

oa _ j:s^sil + 3s)il + zy+^^ 
"^^^^ a? Ts^sil + zf^^^ ' 

As is the so-called Dyer-Roeder parameter that expresses the degree of homogeneity of the 

universe, i.e. As = 1 corresponds to a perfectly homogeneous matter distribution, = 

corresponds to the case in which all the matter is in bound systems. In general, eq. |^ 

can be solved only numerically once the appropriate boundary conditions are specified; in 

a few limiting cases, though, analytical solutions can be found. For the case of interest 

here, in which the matter is supposed to be clumped in a discrete number of collapsed 

objects projected onto the various planes, the relevant limiting case is = for which the 

following solution holds: 

D{zi, Z2) = (1 + zi)[D{z2) - D{zi)l (10) 

where 

D{z)=^r , • (11) 

^0 Ji y\^JW^^ 



This is the formula we use to calculate distances in eq.^ We note that the relation between 
the Dyer-Roeder angular distance and the more familiar luminosity distance, is always 
equal to D = {1 + zy^Di- 

The deflection angle caused by the n-th plane, az„ {Cn) is defined as the gradient of the 
deflection potential entering eq. ^. Throughout this paper we will assume for simplicity 
and for computational economy that lenses are point-like masses, an approximation whose 
limitations are discussed in § 7. The mass distribution of this lens ensamble is derived via 
the PS formalism (eq. ^). It is easy to verify that in this case 

a^Ain) = v^^^ja) = — — i^' (12) 

where is a reference mass and j runs on all the objects on the n-th plane; ^j"'' is the 
position of the j-th mass on the n-th plane, M-""*. In GL studies it is very common to use 
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angular coordinates rather than hnear ones in order to define the deflection of hght rays. 
We follow this tradition by introducing the angular impact parameter 9^ = C,n/Dn, where 
Dn is the distance from the observer to the n-th plane. Eq. ^ then becomes (also using 



eq. [T|) 



where the auxiliary variable // = {c/ Hq){DI^ — Dj^) has been used. From eq. |1^ one can 
appreciate the existence of a characteristic redshift-independent value of 9 which will be 
used as the natural unit angle for the numerical simulations 



^ . 2 X 10-y|.a..csec. (14) 

Because of image splitting and distortion, GL produces a magnification of the fiux from 
a given source. This can be understood by recalling that the angular size of the source 
is increased and, since the surface brightness is conserved, the net fiux received is olso 
increased. In principle, we can calculate the total magnification fi {i.e. including all the 
possible multiple images) of a source located, for example, on the plane 9n iteratively by 
using eq. IT3|: 



fx-\9i) = det 



d9 



N 



(15) 



d9i 

and summing over all multiple images. Regions in the source plane for which fi —>■ oo are 
known as caustics. In practice, this method can be applied only to a very limited number of 
cases in which the lens equation can be inverted, that is 9i can be expressed as a function 
of ^AT. In the next Section we will describe how fi is calculated in practice within the 
numerical scheme used to solve the problem. 
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4. Numerical Simulations: Method and Description 



The problem in which we are interested here, i.e. the magnification of distant sources 
by the intervening cosmological distribution, and outlined in the previous Sections, cannot 
be solved by any analytical technique. Thus, we have to resort to numerical simulations; in 
particular we adopt the so-called ray-shooting method. A full description of this scheme 
can be found, for example, in Wambsganss (1990) (see also Schneider & Weiss 1988); in the 
following we outline only its main features. 

Let us consider a given cosmic solid angle, uj, and study the propagation of light rays 
through it; the main idea of ray-shooting is to follow such propagation in a direction that 
is opposite to the physical one, namely from the observer to the source. This propagation 
from a plane to the next one is mapped by eq. |TB| above for each ray. Thus, illuminating 
all the images produced by the mass distribution contained in tu, we can obtain an accurate 
description of the spatial pattern of the caustics and calculate the source magnification. 

We now determine the cosmological mass distribution to be projected onto the various 
planes. We start by calculating the total mass contained in the cone defined by the solid 
angle uj. The comoving volume of such cone between the two redshifts zi and Z2 is given by 
(hinder, 1988a; Carrol et al. 1992) 



where H{z) is the Hubble constant at redshift z. The corresponding total mass is then 
simply = (^M + ^v) PcV iyJ ^ Z2) , where the critical density pc = {3Hq/8ttG). We 

then subdivide this mass among various point lenses distributed according to the prescribed 
mass distribution function obtained from PS (eq. ^) and appropriate for the cosmological 
model (SCDM, LCDM, CHDM) under study. Specifically, we extract via a Monte Carlo 
procedure values of the masses from such distribution, until their sum exceeds the value 




1 2 



(16) 
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of M-i.j calculated for the two redsliifts bracketing the redshift of the relevant plane. This 
coupled PS+Monte Carlo procedure, being very flexible and handy, can be used to extract 
from any hierarchical cosmological model the information required for the setup of a GL 
numerical experiment. Since the PS formalism does not provide us with any information 
about the spatial distribution of the collapsed structures (we neglect in this paper any GL 
effect related to the possible different distribution of baryons), it is necessary to make an 
ad hoc hypothesis concerning it. The simplest assumption is that the lenses are spatially 
uncorrelated and randomly distributed on the planes. Of course, this neglects the fact that 
in reality a certain, albeit not completely understood, degree of clustering could be present 
even at very high redshift. The limitation of this and the point lens hypotheses, that are 
the main ones in this paper, will be further discussed in §7. 

With this background, we are now ready to describe the main features of the numerical 
computations and of the determination of the magnification probability of a lensed source. 
On the first plane, we define an orthogonal grid made of A*"; x A*"; pixels of size A; and we 
shoot a light ray at each cell center. We then apply the lens equation |13| to the angular 
impact position 6i of each ray to obtain all the subsequent positons 02,6n,(^N up to the 
source plane. There we collect the final impact position in another grid made of A^^ x Ng 
cells of size A^. The choice of A^ must be a compromise to minimize the effects of two 
statistical biases: on the one hand the number of rays collected in each cell has to be large 
enough to allow for a good statistical significance of averages; on the other hand, cells that 
are too large to prevent good spatial resolution of the caustic patterns. For this reason we 
constrain A^ to be 

A; < A3 < L„ (17) 

where Lg = NgAg. In general, Lg < Li = NiAi in order to avoid well-known border effects, 
namely the flux attenuation in cells close to the grid borders due to the neglect of masses 
external to the cone. The final product of a numerical run is a A^s x A^^ square matrix; 
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each element of the matrix represents the number of counts, Uij, in the cell The 
magnification matrix can then be obtained by evaluating the change in the actual number 
of rays per cell in the two grids: 

f^iji^) = ni.iz)-^. (18) 

s 

Due to the statistical character of the lens distribution, the information contained in 
fiij{z) has to be interpreted statistically, as well. This can be done by introducing the 
magnification probability function 

I Ns 

2;) = — ^ 7^ [^ij{z) - ^] , (19) 

^^s i,j=l 

where 7i = 1 if /ijj = fi, and TC = otherwise. An additional useful quantity is the 
cumulative magnification probability function 

p{> /X, z) = E - /^), (20) 

^^s i,j=l 

where Oh is the Heaviside function. Clearly, since fi is calculated for each cell, this implicitly 
implies that the probability distribution is appropriate to point sources. 

We conclude this Section giving the numerical values of the above parameters used in 
the simulations. We use = 50 planes, equally spaced by a redshift interval Az = 0.2, thus 
spanning the total redshift range 2 = — 10. The number of rays followed in each simulation 
is = 1.85 X 10''; they are uniformly distributed over a solid angle u = 2.8 x 10"^ sr, 
corresponding to a 420 x 420 arcsec field. We assume M* = lO^^M© (corresponding to 
the largest mass predicted by PS in each run) and consequently 9^ = 42 arcsec (see eq. 
T¥). Then it is Lg = 76^ = 294 arcsec and, as stated above, Li = 106^ = 420 arcsec; 



As = Lg/Ns = 0.96 arcsec for Ng = 300. Since we are shooting one ray per cell on the 
first plane, the minumum angular resolution is Li/Ni = 0.01 arcsec, a typical deflection 
caused by a mass M ^ 1O^°M0. Mainly for this reason we have applied a lower mass cutoff 
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to the mass distribution = 5 x IO^^Mq which represents a good compromise when 
computational economy is taken into account. Of course, the number of lenses varies with 
redshift; the largest value in a plane is ~ 600 lenses with little difference among different 
cosmological models. Finally, we point out that, although the solid angle uj = 2.8 x 10~^ sr 
considered in the numerical runs is relatively small, the statistical significance of the results 
is guaranteed by the sufficiently large number of masses typically contained on each plane, 
of order of several hundreds. 

5. Results 

Our results are essentially contained in the magnification maps obtained from the 
numerical simulations. In the following, we first present and describe them in general terms, 
underlying the differences among cosmological models; next, we analyze them in detail by 
means of the magnification distribution function described above. 

5.1. Magnification Maps 

Fig. shows three magnification maps relative to the three cosmological models under 
study. Typically, = 50 maps are obtained by gradually increasing the redshift of the 
lensed point source mimicking a distant SN, possibly occurring in a Pop III object. Here 
we show as an example those corresponding to z = 6 for each model, which nevertheless 
are representative of the general lensing features of the various cosmologies. These maps 
give the magnification of the flux source as a function of its putative position in the source 
plane. Note that, since we have normalized all fluctuation spectra such to reproduce the 
observed number of clusters at z = 0, the three models are very similar at low redshift and 
start to differentate as the source is moved back in time. Hence, the magnification patterns 
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at sufficiently high z truly evidentiate intrinsic differences in the mass distribution function 
predicted by the various models, allowing for a meaningful comparison. 

Starting from the map corresponding to the SCDM model, we note that caustics 
appear rather intense and concentrated in structures with sizes of the order of ~ 12 cells, or 
^11 arcsec. About 10 high magnification regions with yU ^ 50 can be identified, which are 
produced by the most massive lenses encountered by the light rays along their cosmological 
path. The median magnification calculated on the entire map is /i = 1.5. The LCDM 
magnification map in Fig.pD, to a ffist sight looks similar to the previous one, again showing 
the same filamentary structures. Intense caustics are less frequent and of smaller size; a 
very strong magnification /i ^ 120 is found at one location, but the median magnification 
value is slightly lower than for SCDM, /2 = 1.47. CHDM models, are instead peculiar. In 
fact, they produce a "grainy" pattern distribution, constituted by a large number of low /i 
caustics of small angular size, of the order of a few grid points (~ few arcsec). Caustics 
more intense than ^ 30 are absent, and the median magnification is /i = 1.4. 

These differences (or resemblances) among various models can be understood easily on 
a qualitative basis, noting that CHDM models predict less massive objects than SCDM or 
LCDM at any z ^ 0.5. This emerges in the different effects of GL which, in the ffist case 
is dominated by the effects of masses typical of galaxies or small groups (M ^ 10^"^ Mq), 
whereas for cold dark matter only models is influenced by the presence of cluster-type 
lenses (M ^ IO^^Mq). The less sensible differences between SCDM and LCDM models 
arise mainly from the fact that for a given redshift the typical lens mass scale in the LCDM 
is about one order of magnitude smaller (see Fig.|l|), whereas the lens number density is 
roughly the same in the two models. This implies lower magnifications and caustic sizes, 
but apparently similar magnification patterns. 
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5.2. Magnification Probability 

Fig. ^ shows the cumulative distribution function z) for the three models. Such 
distributions present a moderate degree of evolution up to ^ ~ 5 (CHDM) and z ^ 7 
(SCDM/LCDM). This behavior, as already stressed before, reflects the different trends of 
the mass function with redshift. All models predict that statistically large magnifications, 
/i ^ 20 are achievable, with a probability of the order of a fraction of percent, the SCDM 
model being the most efficient magnifier. Also the shape of the various curves is rather 
similar. In order to appreciate better the differences among the models and to quantify 
the magnification probabilities in more detail, we introduce two parameters derived from 
the parent distribution P{fi,z). These are: (i) fiio{z), or the value of fi with a probability 
P(yU, z) = 10%, which gives an idea of the magnification of the sources that are more likely 
to be detected, and (ii) P{fi > 10, z), or the probability associated with high magnifications 
(we have set the threshold at /i = 10 to parallel the previous analysis in terms of the 
two-point correlation function), which estimates how likely the detection of highly magnified 
sources will be. The evolution of these two parameters is illustrated in Figs. ^-P. From 
the behavior of fiio^z), we see once again that CHDM models produce (statistically) lower 
magnifications than the other two models. Nevertheless, all cosmologies predict that above 
z ^ A there is a 10% chance to get magnifications larger than 3, which, observationally, 
corresponds to a magnitude gain of Am ^ 1.2. Also, pushing the realm of observations 
to very high z potentially allow discrimination of the models, or at least to discard some 
of them, provided that enough lensed point sources can be detected to reproduce P(/i,z). 
These (and other) observational implications are discussed in the next Section. Similar 
considerations hold also for the evolution of P(/i > 10, z), which shows a steep increase 
of the probability for all models, followed by an almost constant plateau at about 1% 
probability. Thus, in a field of about 25 square arcminutes we expect that statistically, an 
area of at least 30" x 30" will produce magnifications larger than 10 times, corresponding to 
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a magnitude gain of Am — 2.5. Finally, it is interesting to note that at low redshift z ^ 3 
there is a close resemblance in the magnification evolution of LCDM and CHDM models, 
which differentiate only at earlier epochs. 

6. Observational Implications for NGST 

According to hierarchical models of structure formation the first luminous objects in 
the universe should have total masses of order IO^Mq and form at redshift z 20 — 30; 
larger objects then form by merging of these building blocks. However, a time gap between 
the first and subsequent generation of structures could be present, due to the possible 
feedbacks (H2 photodissociation, reheating of the IGM due to stellar energy injection: see 
Shapiro et al. 1994, Ferrara 1998). Thus, for quite a long cosmic time interval these small 
objects could be the only visible trace of galaxy formation at the redshifts that we hope 
to be able to unveil in the near future. Given their small mass, these objects are likely to 
be faint: for a reasonable mass-to-light ratio for young galaxies ~ 0.1, their luminosity 
is 2 X 10^° erg s~^. A Type II SN is typically one hundred times brighter. Thus, for 
periods even longer than a year (taking into account the time stretching oc {1 + z) of the 
SN light curve) the SN outshines its host galaxy. Since Type II SNe originate from massive 
stars, a necessary requirement is that the (unknown) IMF of these objects is flat enough to 
extend into this regime. Clearly, this condition is satisfied in the local universe, but it is 
not necessarily so in the conditions prevailing when the universe was young. Nevertheless, 
it is intriguing to speculate about the observational perspectives to detect very high z SNe, 
a possibility to which our hopes to investigate directly the primeval star/galaxy formation 
are closely tied. 

The advances in technology are making available a new generation of instruments, some 
already at work and some in an advanced design phase, which will dramatically increase 
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our observational capabilities. As representative of such class, we will focus on a particular 
instrument, namely the Next Generation Space Telescope (NGST). In the following we 
will try to quantify the expectations for the detection of high z SNe and the role that the 
gravitational lensing can play in such a search. 

We will assume that NGST (i) is optimized to detect radiation in the wavelength range 
from \min = 1/im to Xmax = 5/im {i.e. J-M bands), and (ii) can observe to a limiting flux of 
^NGST = 10 nJy in 10^ s in that range, which should allow for low-resolution spectroscopic 
follow-up. We also assume that a Type II SN has a black-body spectrum (Kirshner 1990) 
with temperature Tsn and we fix its luminosity Lsn = 3 x 10"^^ (Woosley & Weaver 1986; 
Patat et al. 1994). This constant luminosity plateau lasts for about ^80(1 + 2;) days, after 
which it fades away. 

Fig. 1^ shows the AB apparent magnitude of a SN as a function of its explosion redshift 
in four wavelength bands (J, K, L, M) in the assumed NGST sensitivity range, and we 
compare it with the instrument flux limit. Note that we have taken into account absorption 
by the IGM at frequencies higher than the hydrogen Lya line. The plot also assumes that 
Tsn = 25000 K, a value corresponding to the temperature approximately appropriate to 
the first 15(1 + z) days after the explosion (Woosley & Weaver 1986). After that, the 
temperature drops to T ^ 6000 K, where hydrogen recombines, for the entire duration of 
the plateau (see below). From an inspection of Fig. ^ we deduce that, without considering 
GL magnification effects, the J band offers the best perspectives to observe such an object: 
at the assumed sensitivity limit, NGST should be able to reach z ^ 4. Including GL 
magnification enhances dramatically the observational capabilities. In fact, as one can see 
from the solid line set of curves in Fig. ^ even allowing for a magnification /i = 3 only, 
which in all models has a magnification probability larger than 10% at high z (see Fig. |^), 
this pushes the maximum redshift at which SNe can be detected up to 2; ~ 9 in the K band. 



-22- 



A different way to appreciate the GL effects is illustrated by Fig. ^ There we allow 
for the SN temperature to vary in the large interval Tsn ^ 4000 — 10^ K and we ask which 
is the region of the {1 + z) — Tsn parameter space in which NGST could possibly detect 
SNe in primordial objects. This region is constrained by the requirement that the radiation 
from the SN (or, indicatively, the maximum of its black-body spectrum) falls in between 
Xmin and Xmax witli a flux larger than J^ngst- In the absence of magnification, the upper 
redshift boundary of the detection area (marked as "observable" in the Figure) is 4 < 2; < 8 
for 10^ K < ^ 3 X 10^ K for SCDM and CHDM models (both having = 0) and 
3 < 2 < 5.5 for 8000 K<Tsn^ 2.5 x 10^ K for LCDM. This difference is due to the 
larger luminosity distance in LCDM models. We can now compare this result with the one 
obtained considering the GL magnification of the SN, as calculated in the previous Sections. 
As an indication, we use for /i the value of fiw{z) shown in Fig. |^: this quantity should 
represent a good compromise between the two, usually conflicting, requirements that a high 
magnification and high probability event occurs. In this case, the upper redshift boundaries 
(the three uppermost declining curves in Fig. ^) of the detection area are shifted towards 
higher redshift, in a way that depends on the cosmological models, which are now clearly 
differentiated because of the concomitant effect of their different fiwi^) and luminosity 
distances. Also, the total area is generally increased, thus considerably enhancing the 
detection chances. The maximum redshift at which a SN can hopefully be observed is now 
above z = 10, the limit of our GL simulations, for all models, and it can be extrapolated 
to z ~ 12 for SCDM. The most favorable SN temperature to maximize the depth of the 
search is Tsn ~ 14000 K, but a large range, 10^ K < < 3 X 10^ K, allows to explore 
the universe illuminated by SNe at redshifts larger than eight. This range is limited at high 
Tsn by integalactic absorption of SN photons above the hydrogen Lya frequency. 

So far, we have emphasized how GL magnification can help us detecting very high-z 
SNe. However, magnification maps can also contain precious information to discriminate 
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among different cosmological models. There are at least two ways in which one can think of 
exploiting such information. In principle, the interpretation of SN number counts based on 
the results shown in Fig. ^ could already provide a first insight on the different properties 
of the models, since the 10% magnification probability contours are located at different 
positions. It follows that one would expect different decreasing trends in the number counts 
with redshift above the fi = 1 line. However, this method crucially depends on the star 
formation/SN production history in the universe; hence a firm interpretation of the results 
would require a good knowledge of such quantity. 

A different, more promising, strategy is instead the reconstruction of the (cumulative) 
probability distribution, P{fi,z), shown in Fig. ^. Suppose that, by monitoring a given 
field (the planned NGST field of view has a size 240 x 240 arcsec, very similar to the 
294 X 294 arcsec maps shown in Fig. |^) at some appropriate time interval that depends 
on redshift, a statistically significant number of SNe are discovered following their light 
curves. From the spectral energy distribution, say in the J K L M bands, and assuming 
a black body spectrum for the SN, one can determine the ratio (1 + z)/T. If follow-up 
spectroscopy is available, the redshift (and, consequently, the temperature) can be 
determined straightforwardly. However, even if this measure is not possible, one can try a 
different strategy based on the fact that the light curve is stretched by a factor {1 + z). It is 
well known that Type II SNe fall essentially in two categories, characterized by the presence 
of a luminosity plateau (SNell-P) or by an almost linear luminosity decrease (SNell-L). 
The duration of the plateau in SNell-P has a distribution which is peaked around the mean 
value 70 days. Hence, lacking a better estimate, the redshift can be determined directly by 
properly sampling the light curve. Obviously, the accuracy of this method cannot compete 
with a spectroscopic determination. Nevertheless, there is a fortunate coincidence inherent 
to the properties of the magnification pattern that makes this error not so crucial when 
comparing different cosmological models. In fact, we see that for redshift ,2^5 the P(/i, z) 
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curves approach a "limiting" curve (see Fig. |^), i.e. there is no magnification evolution, 
appropriate for each model. This is due to the fact that in hierarchical models objects 
present at such high z essentially do not contribute to further magnification because of 
their low mass. It follows that, even if the redshift is not very precisely known we can 
still reconstruct P(yU, z) to a very good accuracy by binning together all the detected SNe 
with z ^ 5. This is clearly an advantage of the magnification method with respect to the 
standard direct methods based on the luminosity distance for which a precise determination 
of the redshift is crucial, as far as discriminating among cosmological models is concerned. 
The final step is represented by the determination of the absolute luminosity of the plateau 
phase from which the magnification follows directly. Again we are facing the fact that Type 
II SNe are not proper standard candles: althoug the plateau luminosity is observed to vary 
only by a factor ^ 2, there are exceptions to this rule (Patat et al. 1994). With this note of 
caution, we can have at least a first estimate of the magnification. 

We conclude this Section by mentioning that Type II SNe have also been used as 
distance indicators via the so-called Expanding Photosphere Method (Kirshner & Kwan 
1974, Montes & Wagoner 1995), which has been rather successful in distance determinations 
of closer SNe; this technique has recently been adapted for application to moderate redshift 
SNe (Schmidt et al. 1994a). In short, the luminosity distance (and hence constraints on 
cosmological models) can be obtained from the formula 



Driz) 



where T is the SN's observed color temperature, fi{z)fx is the observed magnified fiux 
density, Bx{T) is the Planck function and Ca is a distance correction factor derived from 
models to account for the dilution effects of scattering atmospheres; the photosperic radius, 
i?, can be obtained from optically thin lines and considering free expansion of the ejecta; 
for this reason, spectral data are necessary (Schmidt et al. 1994b). Thus, this method is 
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suitable for redshift low enough that fj,{z) 1; for the redshifts of interest here, a model for 
the flux magnification has to be nevertheless assumed. 

7. Discussion and Summary 

In this paper we have explicitly calculated the flux magnification of distant sources 
due to gravitational lensing effects of the intervening matter distribution in the universe for 
different cosmological models. We have particularly concentrated on Type II SNe, since at 
very high z they are likely to offer the best detection perspectives, thus giving us a clue on 
the star/galaxy formation history in the young universe. 

We now discuss briefly the approximations and assumptions made. We have assumed 
that the lenses are point masses. Although this approximation is commonly used in 
GL problems and numerical simulations (Schneider & Weiss 1988) due essentially to its 
simplicity, it remains an idealization of the true lens matter distribution. In reality, even 
maintaining the less restrictive spherical symmetry hypothesis, lenses arc more likely to 
be extended objects with a given radial density profile. Often, this is taken to be an 
isothermal profile, since this distribution tends to reproduce reasonably well the properties 
of galaxy clusters (Narayan & Bartelmann 1997) or leads to a good agreement between 
COBE data and CDM cosmological models (Kochanek 1995). It is well known (Vietri 1985) 
that an isothermal sphere has generally a lower magnification power at high magnifications 
than a point-like lens of the same total mass. Thus, the magnifications we find could be 
overestimated by some factor, which depends on the redshift of the source. The magnitude 
of such effect is hard to estimate since we are considering a thick gravitational lens and it 
is not obvious how the magnification is modified with respect to the single object case, for 
which this factor is typically < 10, depending also on the compactness of the lens mass 
distribution. It has to be pointed out that, in any case, a considerable debate is present 
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in the literature on the detailed shape of dark matter halos. CDM models predict that 
halo density profiles are essentially self-similar with a weak mass dependence (Lacey & 
Cole 1994; Navarro et al. 1997); specifically the dark matter density should decrease with 
radius as Ph^^{r) oc at small radii, down to the resolution of the simulation. Recent 
observational work, however, appears to disagree with this prediction. Salucci & Persic 
(1997) have suggested that dark matter halos have a core, i.e. a central region of almost 
constant density, whose size increases with galaxy luminosity both in absolute units and as 
a fraction of the optical radius. 

The second approximation is the neglect of a possible clustering of the lenses. Since 
this information cannot be derived directly from the PS formalism used here, an additional 
hypothesis is necessary to assign the spatial distribution of the lenses on the various planes; 
as a first step we have adopted a random distribution. Again, reality might be different. 
Recent works (Hamilton et al. 1991, Matarrese et al. 1997) have been rather successful 
in reproducing observational data by detailed modelling of the clustering evolution of 
various populations of objects (QSOs and galaxies). These results arc dependent on the 
bias evolution which is still only partly understood (Mo & White 1996, Kauffmann et 
al. 1997, Catelan et al. 1997). A detail study of the effects of clustering of galaxies on the 
magnification distribution function has been carried out by Jaroszyhski (1991). Prom Pig. 4 
of that paper, and from the related discussion, it is seen that the differences with respect to 
a randomly distributed galaxy/lens population are negligible up to values of ~ 30. In the 
high magnification limit, a random distribution leads to a relatively lower probability: for 
example, for ^ 100, P{^) is about 3 times smaller in that case. In a subsequent related 
and more refined work (particularly with a better statistics), Jaroszyhski (1992) instead 
found that the discrepancy between random and correlated lens models are even weaker, 
and their associated magnification probability distributions would be undistinguishable for 
flat cosmological models. The likelihood of observing distant supernovae depends only on 
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the magnification probability. As our results have been calculated assuming a magnification 
with a probability of 10%, typically corresponding to ;U = 2 — 4 (see Fig. it follows 
from the above discussion, that they should not be sensitive to the particular lens spatial 
distribution adopted. Instead, -P(yu) depends more strongly upon the lens mass distribution, 
as it is clear from the differences among the three cosmological models studied here. It is 
interesting to note that in any case as lens correlations tend to enhance high magnification 
probabilities (Jaroszyhski 1992), it will counterbalace to a certain extent the effect of the 
point-like lens approximation discussed above. 

Finally, our description contains two additional hypotheses that should have only 
a minor effect on the results: (i) the baryonic mass fraction is the same in all lensing 
objects (ii) we neglect the possible occurrence of microlensing due to stellar mass objects in 
intervening galaxies. 

To conclude the discussion, we comment on our choice of Type II SNe as the primary 
focus of our efforts. As we have already stated above these events could represent the most 
distant observable ones in the early universe. This statement is based on the recognition 
that a SN is expected to be at least 100 times more luminous than their host Pop III 
object. This is even more true for Type la SNe, that are on average brighter by about 1.5 
mag; moreover. Type la SNe are known to be very good standard candles and, for this 
reason, they are widely used to determine the geometry of the universe (Garnavich 1997; 
Perlmutter 1998). Hence, they could be excellent candidates for the present purposes. 
Unfortunately, it is reasonable to expect that this events at 2; ^ 10 are very rare, since they 
arise from the explosion of C-0 white dwarfs triggered by accretion from a companion; 
this requires evolutionary timescales comparable with the Hubble time at that redshift. A 
different possibility would be represented by Type lb SNe, which originate from short-lived 
progenitors: however, they share problems similar to Type II SNe, i.e. they are poorer 
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standard candles and are fainter than Type la SNe. In any case, our results are valid for all 
three types of SNe since they give the relative flux enhancement due to GL magnification. 

The summary of the main results of this paper, under the assumptions made and 
discussed above, is the following: 

• The magnification patterns are different for the three cosmological models considered: 
for a given z, caustics are more intense and concentrated in SCDM models and their pattern 
becomes more "grainy" for LCDM and, particularly, for CHDM. Median magnifications are 
~ 1.5 for all models. 

• The magnification probability function presents a moderate degree of evolution 
up to z 5 (CHDM) and 2; 7 (SCDM/LCDM). All models predict that statistically 
large magnifications, // ^ 20 are achievable, with a probability of the order of a fraction of 
percent, the SCDM model being the most efficient magnifier. All cosmologies predict that 
above 2; ~ 4 there is a 10% chance to get magnifications larger than 3. 

• We have explored the observational implications of the above results taking NCST, 
with the planned instrumental characteristics, as a reference instrument. We find that 
NCST should be able to detect and confirm spectroscopically Type II SNe up to a redshift 
of 2; 4 in the J band (for Tsn = 25000 K); this limit could be increased up to 2; 9 
in the K band, allowing for a relatively moderate magnification. The detection area in 
the redshift-SN temperature plane depends on cosmological models. When magnification 
is taken into account, and for the most suitable values of Tsn ~ 14000 K, the maximum 
possible detection limit is pushed well above z = 10. Possibly promising strategies to 
discriminate among cosmological models using their GL magnification predictions and very 
high- 2; SNe are sketched. 
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Fig. 1. — Fraction of collapsed objects per unit logarithmic mass for various cosmological 

models with parameters as given in the text: SCDM (solid line), LCDM (dashed) and CHDM 
(dotted) and for two different redshifts: z — 3 (light curves) and z = 5 (thick). 
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Fig. 2. — Schematic view of the Thick Gravitational Lens scheme with symbols defined in 
the text. 




Fig. 3. — Magnification maps for the three cosmological models: (a) SCDM (b) LCDM (c) 
CHDM at redshift z — 6 with symbols defined in the text. 




Fig. 3b - 
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Fig. 4. — Cumulative magnification probability for four different redshifts z — 3 (solid line), 

z = 5 (dotted), z = 7 (dashed), and z = 9 (dot-dashed) for (a) SCDM, (b) LCDM, (c) 
CHDM models, respectively. 
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Fig. 5. — Values of magnification with associated probability larger than 10% as a function 

of source redshift in the three cosmological models: SCDM (solid line), LCDM (dashed), 
CHDM (dotted). 



-43- 




1 23456789 10 11 

1+Z 

Fig. 6. — Probability associated with magnifications higher than /i = 10 as a function 

of source redshift in the three cosmological models: SCDM (solid line), LCDM (dashed), 
CHDM (dotted). 
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Fig. 7. — AB apparent magnitude of a SN as a function of its explosion redshift in the four 
wavelength bands J, K, L, M (from the uppermost to the lowermost curve) for the SCDM 
model. The case for a moderate (/i = 3) magnification (solid curves) is compared to the 
one in which magnification is neglected (dotted). The NGST flux limit is also shown; the 
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Fig. 8. — SN detection area (marked as "observable") in the redshift-SN temperature, Tsn-, 
plane neglecting GL magnification (curves /i = 1) and with magnification /iio(2) for the 
three cosmological models: SDCM (sohd curve), LCDM (dashed), CHDM (dotted). Note 
that for /i = 1 the SCDM and CHDM curves overlap. The area is also bound by the lower 



